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Abstract 

This paper considers the preliminary development of a general optimization pro- 
cedure for tetrahedron formation control. The maneuvers are assumed to be impulsive 
and a multi-stage optimization method is employed. The stages include (1) targeting 
to a fixed tetrahedron location and orientation, and (2) rotating and translating the 
tetrahedron. The number of impulsive maneuvers can also be varied. As the impulse 
locations and times change, new arcs are computed using a differential corrections 
scheme that varies the impulse magnitudes and directions. The result is a continuous 
trajectory with velocity discontinuities. The velocity discontinuities are then used to 
formulate the cost function. Direct optimization techniques are employed. The pro- 
cedure is applied to the NASA Goddard Magnetospheric Multi-Scale (MMS) mission 

to compute preliminary formation control fuel requirements. 
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Introduction 


The purpose of this paper is to introduce and to motivate further study of some 
of the ideas and concepts needed to successfully fly the NASA Goddard Magneto- 
spheric Multi-Scale mission. Specifically, this paper reports the results of a prelimi- 
nary study on the development of a general optimization procedure for tetrahedron 
formation control. In general, for this type of in-situ mission, maintaining a tight 
formation during the orbital evolution is not necessary. It is of importance, however, 
to be able to change the tetrahedron size and to maintain shape metrics (quality 
factors [1]) within “acceptable” values during specific times. Typically, the quality 
factor maximization can be achieved by a judicious choice of initial conditions. The 
formation then evolves naturally (without maneuvers) while collecting science data. 
In this investigation, the maneuvers needed to transfer the formation from its current 
configuration to the initial conditions at the beginning of the science arc are exam- 
ined. In an effort to minimize fuel expenditure, a multi-stage optimization method is 
employed. The algorithms developed for this investigation do not make any assump- 
tions on the location and number of burns that will be used. Direct optimization 
techniques that have been presented and examined previously are used [2j. No claims 
are made about the global optimality of the solutions, only that the solutions com- 
puted are local optima given the inputs and constraints. In fact, an experienced 
analyst might be able to take some of the resulting solutions and further improve 
them via numerical experimentation. Nevertheless, the ultimate goal is to automate 
the mission design process as much as possible by quickly computing solutions that 
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are at least local minima. The methodology is applied to the MMS mission to gain 
understanding about the fuel needed to perform certain tasks. 

Previous Work 

Tetrahedron formations are excellent scientific platforms for magnetic field stud- 
ies. A tetrahedron mission, Cluster II, is currently flying and operating success- 
fully [3]. In the process of planning the Cluster mission, fuel optimization of the 
maneuvers needed to initialize, modify and maintain the formation was considered 
by Rodriguez-Canabal and Bello-Mora [4]. Later a different optimization method was 
employed by Schoenmaekers [5]. In Schoenmaekers’s paper, both methods are com- 
pared for one mission scenario. Moreover, a further-developed strategy is presented 
by Bello-Mora and Rodriguez-Canabal [6]. Implementation and operational results 
for the maneuvers are reported by Hockens and Schoenmaekers [7]. These Cluster 
references provide a solid foundation for the work in the current investigation. 

In this paper the mission of interest is the MMS Mission. MMS will determine the 
small-scale basic plasma processes which transport, accelerate and energize plasmas in 
thin boundary and current layers. These processes control the structure and dynamics 
of the Earth’s magnetosphere [8]. For the interested reader, preliminary mission 
design and analysis for MMS, which includes a double lunar swingby for one of its 
phases, has been presented by: Edery and Schiff [9], Edery [10], Guzman and Edery 
[11,12], Petruzzo [13], and Hughes [14,15]. 
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Approach 


For the purpose of this investigation the formation flying design is divided into 
two distinct arcs. The first arc includes a sequence of maneuvers that transfers the 
formation from its current configuration (position and velocity states) to a configura- 
tion appropriate for the start of the second arc, the science one. Denote the formation 
configuration at the beginning of the science arc as the target configuration. See Fig. 
. Note that the science arc does not contain any maneuvers. In fact, the absence of 
maneuvers benefits the science data gathering by limiting the orbital disturbances. 
The two arcs then lead to two coupled optimization problems: (i) compute the initial 
states of the science arc that maximize a certain tetrahedron quality factor along the 
science arc, and (ii) compute the maneuver sequence that minimizes the fuel needed 
to transfer the formation from its current configuration to the target configuration. 
In this paper, only the second problem is considered. In references [13, 14] the first 
problem is considered. Since the two problems are coupled, the goal is to eventually 
merge these two optimization processes. For now, a target configuration is assumed 
to develop algorithms that solve the second problem. In other words, the four states 
that comprise the target configuration are considered as given for the maneuver op- 
timization process. Nonetheless, as a way of finding the greatest lower bound on fuel 
consumption, rotations and translations of the target tetrahedron have been alfowed. 
Eventually, both quality factor and fuel consumption considerations will determine 
the target configuration. 
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The Concept of a Reference Path 


The trajectory design strategy for MMS involves obtaining a reference or nominal 
path that traverses the magnetospheric regions of interest. Relative to this defined 
reference path, a tetrahedron formation is established at certain times (specified by 
the scientists) during the mission. The science mission consists of four phases: the 
first two phases traverse regions of the magnetosphere close to the magnetic equato- 
rial plane, the third phase contains a double lunar swingby sequence (DLS) with an 
investigation of the deep tail magnetospheric region, and the fourth phase explores 
regions perpendicular to the ecliptic. The different MMS phases and their corre- 
sponding scientific goals are explained by Curtis [8]. The science goals, in turn, help 
to specify the orbital requirements for each phase. The main orbital requirements are 
shown in Table (where R® is the mean Earth radius). It is likely that some flexibility 
in the location of the formation relative to the reference path might be allowed (if the 
magnetospheric regions of interest are still traversed by the formation). Furthermore, 
the reference trajectory could be redefined. Still, since one of the phases includes a 
double lunar swingby, monitoring the formation relative to the reference path is of 
interest. 3 

Coordinate Systems 

For the purpose of this investigation, three orthogonal Cartesian coordinate systems 

are employed. The first system is the geocentric inertial (GCI) frame. The inertial 
3 Because the reference path has been designed with the proper phasing for the lunar encounters. 
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frame is defined using the mean equatorial plane and the mean vernal equinox at Ju- 
lian year 2000. The second frame utilizes the reference path orbital velocity, binormal 
and normal directions (VBN). The third frame is introduced to consider tetrahedron 
translation and rotations. It is placed at the tetrahedron centroid and it is initially 
coincident with the second frame. See Fig. . Table summarizes the frames employed 
(and their respective origins) and introduces the unit vector notation used in Fig. . 

Initial Conditions for the Science Arc 

The algorithms developed for this study are independent of the initial states selected 
for the science arc; yet the actual computations depend on them. Thus, in this section, 
an approximation to select the targets states at the beginning of the science arc is 
explained (other choices are possible). A regular tetrahedron 4 is used to meet the 
science data collection requirements. A reference orbit is specified and the individual 
spacecraft locations (position components) are set using spherical coordinates in the 
VBN frame [17]. 

The velocity directions are kept the same as that of the reference orbit, that is, the 
VBN frame velocity direction. This assumption is adequate when the spacecraft-to- 
reference separation is “small” relative to the reference-to-Earth separation (it should 
be revisited for large separations). Nonetheless, changing the velocity directions can 
be utilized to maximize a certain tetrahedron quality factor. This quality factor opti- 
mization has been investigated by Petruzzo [13]. The investigation has been extended 
4 A regular tetrahedron is useful in cases where sampling of the structure of a field is not as 
important as understanding its transient or fluctuating events [16]. 
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by Hughes [14] to consider both the initial positions (non-regular tetrahedrons) and 
velocities. 

The velocity magnitude is obtained using the two-body problem energy or vis-viva 
equation: v = \J \iiflfr — 1/a), where v is the velocity magnitude, r is the position 
magnitude (radial distance), fi is the gravitational parameter and a is the semi-major 
axis. For orbital periodicity, the reference semi-major axis is utilized for all the 
spacecraft. 

Initial Guess for the Maneuver Sequence 

The first step in optimizing a particular problem is to obtain an initial guess. In this 
case an initial guess is obtained by propagating from some initial configuration. The 
initial states for the maneuver arc are also specified as explained in the Initial Condi- 
tions for the Science Arc section. (This specification is for computational convenience. 
During spacecraft operations the states would be set by the formation configuration 
at the specified start time for the maneuver arc). Then, the reference and the four 
spacecraft formation are propagated without any maneuvers. While propagating, dis- 
crete states along the path are saved. The discretization interval is selected by the 
user. In this case, a discretization in terms of true anomaly (TA) along the refer- 
ence path is utilized. See Fig. for true anomaly discretizations of 180, 90, 45 and 
22.5 degrees. These discretizations provide 3, 5, 9, and 17 “patch states” respectively. 
These patch states are discrete states where maneuvers might be implemented. Other 

anomalies (e.g. mean or eccentric) could be used for the discretization. 5 Now, it is 
5 In fact, the resulting optimal maneuver patterns should be compared in future studies. 
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important to realize that a spacecraft path with, for instance, 17 maneuvers is hard to 
implement in practice (operationally). Nevertheless, numerically implementing such 
a path will show the best locations to perform the maneuvers. Furthermore, many of 
the resulting maneuvers are “small” and can be eliminated and/or combined. As a 
result, the mission analyst can compute operationally feasible solutions. 

Computing the Optimal Maneuver Sequence 

Once the initial and final states have been set for the maneuver sequence, the tar- 
get/optimization process can proceed. The optimization approach involves combining 
the fuel optimization of each spacecraft with the minimization of the total fuel usage 
to achieve a target configuration. The total fuel usage is the sum of the totals for 
each of the four spacecraft (S/C). In terms of the spacecraft engine, only impulsive 
maneuvers are considered (thus only velocity change or AV numbers are presented). 
For computational speed purposes, the optimization is performed sequentially (sepa- 
rately) in two stages: (1) each spacecraft trajectory is optimized by minimizing the 
fuel used by each spacecraft to achieve its target location in the tetrahedron, then, (2) 
the final tetrahedron configuration is varied (translated and rotated relative to some 
arbitrary configuration in the vicinity of the reference path) while minimizing the sum 
of all the spacecraft A Vs. Better" results might be obtained by completely embedding 
step (1) in each iteration of step (2) or by performing step (2) first and then step 
(1). Nevertheless, if the target orientation is specified by the scientific requirements 
and/or by the optimization of a certain quality factor, step (2) might not be allowed. 
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For both stages, as the impulse locations and times change (i.e., as the initial patch 
states in Fig. change), new arcs are computed using a differential corrections scheme 
that varies the impulse magnitudes and directions. See reference [2] for more details. 

In the first optimization stage and for each spacecraft, the transfer trajectory 
problem becomes essentially a rendezvous problem. Specifically, each spacecraft path 
is discretized into a set of n patch states. The independent variables chosen are the 
changes in the internal and final patch state locations and times (i.e., four independent 
variables at each patch state). Thus, there are 4 x (n — 1) independent variables per 
spacecraft. For numerical purposes, the changes are constrained to stay within some 
user defined limits. Moreover, less flexibility is allowed in the final states (position 
and time are constrained to be within 10% of the desired inter-spacecraft distance and 
within 0.5 seconds of the desired final time respectively). Each spacecraft trajectory 
is optimized independently using the total AV as the cost function. The optimization 
method selected is the Sequential Quadratic Programming (SQP) method. This direct 
method is employed due to its efficiency in solving a large variety of cases. (See Ref. [2] 
for a comparison of different methods). 

In the second optimization stage, the target tetrahedron is allowed to translate 
within some limits and to rotate about its initial centroid. The sum of all the space- 
craft A Vs is the cost function. The optimization method selected for the second stage 
is a Genetic Algorithm (GA) that varies 6 independent variables (3 variables for the 
translation and 3 variables for the rotation). Although other choices are possible, this 
stochastic method is selected to avoid local minima and, most importantly, to prepare 
for augmenting the cost function with science and operational constraints (the aug- 
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mentation might result in discontinuities in the cost function). For this preliminary 
investigation, following some empirically developed guidelines [18], a population size 
of N = 46, where 6 is the number of bits in the binary chromosome, and a mutation 
probability of P mut = (6 + \)/2Nb were utilized. These and other details of the setup 
for the GA are still under investigation. 

Some Results for MMS 

The methodology described in this paper is applied to the MMS mission. The 
reference orbital states considered are consistent with the orbit requirements of MMS 
(see Table ). Now, although the algorithms developed for this study are independent 
of the particulars of the force model and integrator, the actual computations depend 
on them. Thus, the dynamical model that is adopted to represent the forces on 
the spacecraft includes the gravitational influences of the Sun, Earth and Moon (all 
obtained from the Jet Propulsion Laboratory Definitive Ephemeris 405 file). Solar 
radiation pressure (flat plate model) and the Earth’s J 2 Earth gravity harmonic are 
also included. For the numerical integration scheme, a Runge-Kutta-Vemer 8(9) 
integrator is utilized. Next, some results and discussion are presented for each mission 
phase. 

Phases 1 and 2 

As a test, different cases (with 3, 5, 9, and 17 maneuvers per spacecraft) are examined 
for one orbit revolution. That is, the initial states are specified at the apogee of the 
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reference orbit and the target states are at the next apogee. 6 Then, for both phases, 
three initial and three final tetrahedron sizes (inter-spacecraft separations between 
each spacecraft) are considered: 10, 1000, and 2000 km. These separations are con- 
sistent with preliminary science specifications [8]. For each initial separation, three 
target separations are considered while keeping the pre-maneuver (initial) conditions 
the same. Therefore, when the initial and target sizes are the same, the problem is a 
formation maintenance problem; otherwise, it is a re-sizing problem. Recall that the 
reference orbital requirements are shown in Table . For now, Phase 2 is assumed to 
have a 10 degree equatorial inclination. 

The results for Phase 1 are in Table and the results for Phase 2 are in Table . 
The majority of the cost function values associated with Phase 1 are higher than for 
Phase 2. Therefore, a more extensive analysis should be performed to understand 
the perturbation effects (e.g. Earth’s J 2 , luni-solar perturbations and solar radiation 
pressure) on the tetrahedron evolution during the different orbital phases. Another 
interesting aspect is that the AV expenditure between the four spacecraft exhibits 
large differences in some of the cases. This particular problem can be remedied by 
adding constraints that require that the total AV of each spacecraft be within some 
(user defined) range relative to the total AV of any other spacecraft. 

While computing some of these test cases, it was observed that the tetrahedron 
translation and rotations did not help to lower the total AV in the one-orbit rev- 
olution maintenance cases. This fact is not surprising given the fact that only one 
6 At the moment, regions near apogee are of interest to the scientists. However, other regions 
could he considered later on. 
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revolution was considered and the perturbation effects did not cause the formation to 
deviate significantly. Nonetheless, during actual mission operations, more than one 
orbit revolution will most likely elapse before tetrahedron maintenance is required. 
The actual maintenance schedule will depend, among other things, on how much 
degradation of the quality factors can be tolerated while still being able to perform 
“effective” science data collection. 

It is also observed that re-sizing cases such as a 10 to 1000 km and 1000 to 10 km 
have similar costs. It is important to understand that these cases are not the same. 
That is, a tetrahedron with inter-spacecraft separation of 10 km and propagated for 
one orbit revolution undergoes different orbital perturbations than one with a 1000 
km separation. Again, a more detailed perturbation analysis should be performed. 
Note that some cases did not “converge”. This label means that either the differential 
corrector could not meet the requested tolerance or that the AV computed exceeded 
some user defined limit. In Tables and the cases that did not “converge” are all 
three-impulse re-sizing cases that exceeded the user defined AV limit. For these 
non-convergent cases, the initial angle between any two impulses is 180 degrees with 
the actual locations being at the apses. As a result, the differential corrector (which 
utilizes the orbit state transition matrix) encounters a singularity [19]. Specifically, 
after some iterations the fourth spacecraft has a significantly different orbit plane 
relative to the other spacecraft (after the second impulse at perigee). Fortunately, 
this particular problem has physical meaning and constraining the transfer plane can 
remedy the situation. Future work should correct these and other non-convergent 
cases and improve the automation process. 
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Phase 3 


In Phase 3, the formation will traverse the deep tail of the magnetosphere during a 
double lunar swingby sequence [11]. At this point, there is no plan to perform any 
deterministic maneuvers for formation flying. That is, the final maneuver at perigee 
before the first lunar swingby produces a trajectory that requires no further maneuvers 
until Phase 4 begins after the second swingby. Thus, the formation will evolve under 
its natural dynamics in the perturbed Sun-Earth-Moon system. Non-deterministic 
maneuvers, however, will be used for any needed corrections. 

Phase 4 

Next, after the double lunar swingby sequence, it is of interest to know how much fuel 
would be required to restore the tetrahedron at the next apogee. In reference [11], 
three maneuvers (after the double lunar swingby) were used to restore the tetrahedron 
at the target apogee: the maneuver at perigee (TA = 0 degrees), an additional ma- 
neuver at the semi-latus rectum radial distance (TA = 90 degrees), and a maneuver 
at the target apogee (TA = 180 degrees). The lowest AV cost computed was for an 
initial separation of 10 km at the beginning of Phase 3 and a target separation of 200 
km at the apogee after the double lunar swingby. The AV results in reference [11] 
were high because each spacecraft was incorrectly constrained to be at its own apogee 
at the end of the sequence. In this investigation, only the reference spacecraft is at its 
apogee at the target time. Also, the maneuver locations and times can now change. 
See Table for the results. It should be remarked that this table is different from the 
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previous tables for Phases 1 and 2. Specifically, the formation was initialized before 
the DLS. Therefore, the start of the maneuver sequence is at the perigee after the 
DLS. At this perigee, the formation is not in a regular tetrahedron but in some elon- 
gated shape with an average inter-spacecraft distance of 24 km. Also, the formation 
traverses only half of the reference orbit. 

Further Fuel Cost Reduction 

Further reduction in the fuel cost might be obtained in some cases by running 
the optimization stages completely nested. That is, for each new target orientation 
compute the optimal trajectories for each spacecraft. This process requires more com- 
putational resources but can be done with selected cases. Other resources include, 
but are not limited to: (a) varying the final (target) time, (b) establishing the tetra- 
hedron relative to one of the four spacecraft (thus, one spacecraft is not required to 
maneuver), and (c) allowing several orbit revolutions before performing any maneu- 
vers (e.g. maintenance). In option (b), the analyst should check that deviations from 
the desired reference path do not have adverse science and trajectory effects. 

Additional insight is obtained by looking at the required configuration changes 
in terms of the orbital elements. This insight might lead the analyst to lower cost 
solutions. In fact, by computing the orbital elements, it can be shown that formation 
re-sizings require mostly line of nodes and line of apsides changes. Therefore, the 
analysis by Lawden is relevant [20], Furthermore, for MMS, it has also been observed 
that Lawden’s primer vector theory (see reference [21] for example) provides better 
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initial guesses for the strategy examined in this paper. As a result, the primer vector 
approach is also currently under investigation [22]. 

Summary 

The purpose of this paper is to introduce and to motivate further study of some of 
the ideas and concepts needed to successfully fly the NASA Goddard Magnetospheric 
Multi-Scale mission. Specifically, the problem of maneuvering a tetrahedron forma- 
tion while minimizing fuel expenditure is considered. The optimization is successfully 
performed sequentially in two stages: (1) each spacecraft trajectory is optimized by 
minimizing the fuel used by each spacecraft to achieve its target location in the tetra- 
hedron, and (2) the final tetrahedron configuration is varied (translated and rotated 
relative to some arbitrary configuration in the vicinity of the reference path) while 

minimizing the sum of all the spacecraft AVs. The methodology is applied to the 

/ 

MMS mission to gain understanding about the fuel needed to perform certain tasks. 
The ultimate goal is to automate the mission design process as much as possible by 
quickly computing solutions that are at least local minima. Future work includes 
the simultaneous optimization of the maneuver and science arcs, the application of 
the methodology to different mission scenarios, the inclusion of operational consid- 
erations/constraints and the integration of the developed software into the NASA 
Goddard mission analysis tools. 
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Figure Captions 


Fig. 1: Maneuver and Science Arcs 
Fig. 2: MMS Frames 

Fig. 3: True Anomaly Discretization: Initial Maneuver Locations 
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Table Captions 


Table 1: MMS Orbital Requirements 
Table 2: Coordinate Frames 

Table 3: AV Cost as a Function of Initial and Final Separations in Phase 1 
Table 4: AV" Cost as a Function of Initial and Final Separations in Phase 2 
Table 5: AV Cost for Restoring the Tetrahedron in Phase 4 
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Mission Phase 

Perigee 

(*») 

Apogee 

W») 

Semi-Major Axis 
(«•) 

Eccentricity 

Inclination 

(degrees) 

Phase 0 

r 1.2 

12 

6.6 


28.5 10 * 

Phase 1 

1.2 

12 

6.6 


10* 

Phase 2 

1.2 

30 

15.6 

wBm 

10-20* 

Phase 3 

NA # 

> 100 

NA 

NA 

NA 

Phase 4 

10 

40 

25 

0.600 

~ 90** 


* Not Applicable, * Equatorial Inclination, ** Ecliptic Inclination 
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T yP e 

Unit Vectors 

Origin 

Inertial 

Rotating 

Rotating 

(* 1 > hi *3) 

(V, B, N) 

(hi hi h) 

© - Earth’s center 
0 - Reference spacecraft (fictitious) 
C - Tetrahedron centroid 
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Impulses per S/C 

AVsci 

AVsC2 

AVscs 

A' 

lO to 10 km 

3 

0.253 

0.235 

0.077 

0 . 

5 

0.028 

0.023 

0.010 

0 . 

9 

0.029 

0.029 

0.031 

0 . 

17 

0.042 | 

0.036 

0.046 

0 . 










38.559 

38.514 


10 to 1000 km 


41.117 

39.042 

42.195 

67.829 


82.772 

69.385 

77.264 

78.616 


35.835 

39.184 


58.301 

52.567 

47.223 

47.910 


66.300 

59.247 

60.079 

39.409 


10 to 2000 km 


116.749 

120.594 


132.839 

93.236 


61.698 148.880 

146.109 144.089 


1000 to 10 km 


56.586 

32.298 

68.365 


126.301 

74.064 

74.454 


1000 to 1000 km 



1.275 

1.231 


0.592 

0.554 

0.563 

0.754 


1000 to 2000 km 




58.098 

68.990 


2000 to 10 km 



2000 to 1000 km 


45.998 

52.814 


0.752 

0.084 

0.115 

0.156 


N/A 

207.442 

181.795 

223.512 




N/A 

409.516 

361.906 

443.269 


58.320 

69.079 

* 

N/A 

52.718 

59.994 

56.589 

207.885 

50.819 

58.251 

44.049 

188.954 

62.433 

69.373 

60.850 

231.840 


2.457 

2.293 

2.299 

3.189 


♦ 

N/A 

57.785 

208.327 

44.388 

188.736 

61.490 | 

231.939 


144.820 

* 

N/A 

121.030 

114.301 

417.326 

118.155 

88.540 

380.153 

139.962 

122.593 

466.950 


71.591 

* 

N/A 

62.622 

57.765 

211.715 

60.888 

44.602 

192.637 

71.842 

61.705 

236.450 
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Impulses per S/C 


AV$C3 


10 to 10 km 


AVsc4 ! AVtotal [m/s] 



0.043 

0 . 


0.024 

0.153 

0.022 

0 . 

1 

0.023 

0.092 

0.027 

0 . 


0.025 

0.105 

0.037 

0 . 

040 

0.037 

0.151 


23.479 

23.479 

* 

N/A 

19.530 

15.701 

20.617 

19.958 

21.954 

8.847 


9.329 

10.069 

13.354 

57.077 


35.130 

45.762 

21.847 

21.790 


10 to 2000 km 


46.833 46.612 * N/A 

44.764 32.428 28.986 151.940 

23.541 52.335 16.014 113.738 

44.021 16.849 31.597 114.257 



17.165 

14.823 

12.349 


1000 to 10 km 


23.630 
19.980 
16.424 
16.599 I 20.483 


1000 to 1000 km 


22.037 

12.669 

12.704 



3 

0.496 

5 

0.453 

9 

0.302 

17 

0.274 



8.440 

5.613 



0.340 0.688 


1000 to 2000 km 


22.728 23.067 

19.066 20.522 

15.889 19.450 

17.150 16.283 


60.936 


0.308 

0.250 

0.216 

0.216 

2.962 

2.124 

1.563 

1.518 


* 

22.293 

12.596 

22.126 

N/A 

77.494 

60.718 

72.534 


* 

44.310 

25.438 

25.498 

N/A 

156.664 

122.848 

122.745 




78.732 

30.032 

24.644 

22.281 


12.192 

11.007 


1.003 

0.867 

0.563 

0.492 


2000 to 10 km 


115.462 | 144.820 
43.167 
40.185 
33.145 I 41.822 


2000 to 1000 km 


42.702 
22.848 
16.666 20.923 

16.854 I 21.711 
2000 to 2000 km 


265 3.331 

803 2.274 

621 1.607 

612 I 1.575 
did not converge 


22.391 

12.902 

12.933 


0.655 

0.482 

0.407 

0.395 


62.682 

62.505 


6.253 

4.426 

3.198 

3.073 
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